Size of the stable population in the Penna 
bit-string model of biological aging 



K.Malarz, M.Sitarz, P.Gronek, and A.Dydejczyk 

Department of Applied Computer Science, Faculty of Physics and Nuclear 
Techniques, AGH University of Science and Technology 
al. Mickiewicza 30, PL-30059 Krakow, Poland 
malarzQagh . edu . pi 



Abstract. In this paper the Penna model is reconsidered. With com- 
puter simulations we check how the control parameters of the model 
influence the size of the stable population. 

1 Introduction 

The Penna model [1] is popular and efficient way to simulate biological aging 
in genetically heterogeneous population (see [2-5] for review). It is based on the 
mutation accumulation theory which claims that random hereditary deleterious 
mutations accumulate over several generations in the species genome [6]. 

The model is governed by a set of parameters such as environmental capacity 
N,nax SO that the actual population N < Nmax- The birth rate B is the number 
of offsprings from a single parent per evolution time step, providing the parent 
reached age R and is still below maximum reproduction age E. Other parameters 
which control the growth of the population are: number of mutations M injected 
into the baby's genome, apart from the already inherited, and threshold value T 
of active bad mutations above which items are eliminated for this only reason. 

In this paper we would like to check how model control parameters (described 
in the next section) influence the size N{t) of the stable population for long times 
t ^ CO when mutational meltdown [7-12] is avoided. 

2 Model 

In the Penna model [1] each individual in the population is characterised by its 
genome — the Nbn-long binary string. Time is measured by a discrete variable 
t. In each time step {t ^ t + 1), for each individual of age a, the number of bits 
set to one in the flrst a positions i in genome is calculated. Bit '1' represents a 
bad mutation activated when its position 1 < i < a, however mutations above 
a < i < Nut are not harmful. If the number of had and active mutations is 
greater or equal to the threshold value T an individual dies due to too many 
deleterious mutations and so called genetic death takes place. 

Individuals competite among themselves for food and territory: each of them 
may be removed from the population with probability N{t)/Nmax, where Nmax 



represents the maximal environmental capacity. In other words, we introduce 
the Verhulst's factor in order to avoid exponential growth of the population size 
to infinity. The number of bits in genome restricts also maximal age of 
individuals. 

We start the simulation with Nini individuals which genomes contain only 
'0' (no bad mutations). The population reproduces asexually. If the individual is 
older than minimum reproduction age R it is able with the probability b to give 
B offsprings. The offspring's genome is a copy of parent's one, exposed during 
replication to M harmful mutations. Each mutation occurs with probability 
m at randomly chosen position in the genome. We also introduce a maximum 
reproduction age E, so the individual older than E do not clone itself any more. 

3 Results 

We would like to refer our results to the standard one (labelled as 'std' in all 
figures) which we assume to be for Nbn = 32, N^ax = 10^, n{Q) — N{Q)/Nmax = 
0.1, i? = 8, = 32, T = 3, M = 1, B = 3, m = 1.0 and 6 = 1.0. The 
influence of the Nut on population characteristics was discussed earlier in [13, 
14]. The maximal environment capacity Nmax is fixed during simulation. The 
initial concentration of individuals in respect to maximal environmental capacity 
do not influence the results (Fig. 1(a)). 

As one may expect the larger probability of reproduction < 6 < 1 gives 
the larger population size (Fig. 1(b)). The population size decreases when the 
probability of mutations < tt?. < 1 increases (Fig. 1(c)). 

The same situation is for integer variables B and M (see Figs. 1(d) and 1(e), 
respectively). If the number of mutations is too large (i.e. M > 7) the mutational 
meltdown occurs. Also too large value of the number of offsprings (e.g. B > 20) 
may lead to catastrophe: too rapid growth of the population may cause that in 
one time step the environmental maximal capacity Nmax will be exceeded. As 
a consequence of that the population will vanish only due to Verhulst's factor 
and not the genetic one. In our implementation products Mm and Bb may be 
considered as the average number of mutations per individual's genome and as 
the average number of offsprings given every year per reproductive individual 
(i.e. in age R < a < E), respectively. 

In Fig. 1(f) the influence of the threshold T on the population size is pre- 
sented. Decreasing the threshold T results in increasing of the stable population 
size. Setting T = Nbn eliminates genetic death and again only Verhulst's factor 
determinates N{t — > oo). 

Delay in starting the reproduction decreases the population size as presented 
in Fig. 1(g). Obviously, setting R > Nbu destroys population because individuals 
never replicate themselves. For R = 1 — when all newly born babies become 
adults as quickly as possible — a very strong oscillations appear (not shown). 

Fig. 1(h) shows the influence of the maximal reproduction age E on the 
stable population size. For E = 8(= R), when individual reproduce only once 






Fig. 1. The influence of the model control parameter (a) n(0), (b) b, (c) m, (d) B, (e) 
M, (f) T, (g) R, and (h) E on the population size. 



(semelparous organisms [23]) oscillations in the population size appears (not 
shown) . 

4 Conclusions 

The population size of the stable population in the Penna model was already in- 
vestigated for a certain purpose, e.g.: showing advantage of sexual reproduction 
on asexual one [15-18], studying sexual fidelity of males vs higher reproduc- 
tion rate [19], modelling fishing and /or hunting [20-22], explaining mystery of 
Pacific salmon senescence [23], analysing pray-predator systems [12], recovering 
demographic data [13, 17, 24], modelling the oldest old effect [25], modelling con- 
sequence of parental care [26, 27], studjdng the individuals fertility [28], studying 
the migration effect [29], and others [30-40]. 

In this paper we confine our attention to the basic features of the model. 
Main results were collected in Figs. 1. We may conclude the following: 

1. Different time scales must be accounted for, if the 'stable' refers to different 
properties observed in the computer experiment; 

— total population N{t) shows values that do not change significantly after 
about one hundred iteration steps, 

— however, we need at least a couple of thousand steps for age distribution 
[14, 25] and, especially, bad mutations distributions in a genome position 
[12, 15], to become stationary. 

2. Obviously for smaller birth rates Bb the population is smaller. However, 
unlike the logistic model for which N{t -> oo) = BbN,„ax/il + Bb) [30, 
40], the Penna model shows a critical Be below which the population is 
extinct. This is due to the genetic death as result of pumping in more and 
more mutations M into each new generation, and so Bc{M) with the limit 
Be = for M — 0, the logistic case. 

3. When the number of bad mutations M injected into baby's genome is scaned 
from M = (the logistic model), the upper limit Mc for the extinction 
strongly depends on fertility: 

— if i? is smaller, Mc is also smaller, and 

— for sufficiently large B, Mc ^ oo as the high birth rate gives a chance 
to new population (however rich in bad mutations) to reproduce before 
it dies out soon on reaching the threshold T. This is so if the minimum 
reproduction age R is sufficiently small to allow for the reproduction. 

4. The upper age limit E for reproduction influence the stable population N 
and larger E makes bigger N, as it may be expected. This time, however, 
we do not observe any critical behaviour when we scan E. 
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